The energy spectrum of cosmic-ray electrons at TeV energies 
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The very large collection area of ground-based 7-ray telescopes gives them a substan- 
tial advantage over balloon/satellite based instruments in the detection of very- high-energy 
(>600 GeV) cosmic-ray electrons. Here we present the electron spectrum derived from 
data taken with the H.E.S.S. system of imaging atmospheric Cherenkov telescopes. In this 
measurement, the first of this type, we are able to extend the measurement of the elec- 
tron spectrum beyond the range accessible to direct measurements. We find evidence for a 
substantial steepening in the energy spectrum above 600 GeV compared to lower energies. 

PACS numbers: 95.85.Ry 



In stark contrast to hadronic cosmic rays (CRs) 
the lifetime and hence propagation distance of CR 
electrons in the very-high-energy regime is severely 
limited by energy losses via synchrotron radiation 
and inverse Compton scattering. The lifetime of a 
very-high-energy electron can be expressed as: t « 
5 X 105(£;/lTeV)-i((B/5/^G)2 -I- 1.6(u;/leVcm-3))-i 

years, where w is the energy density in low frequency 
photons (hiy ^0.1 cV) in the interstellar medium and 
B is the mean interstellar magnetic field. In standard 
diffusion-dominated models of Galactic cosmic-ray trans- 
port this implies that the sources of TeV electrons must 
be local (< 1 kpc distance) as discussed in e.g. Q and 
A second consequence of these energy-dependent 
losses is that the electron spectrum is steeper than that 
of the hadronic CRs (~£'~^-3 of. E^'^-''). AU measure- 
ments so far have utilized balloon or satellite borne in- 
strumentation (see Q for a review). However, the rapidly 
declining flux makes such direct measurements at high 
energies difficult. It has been suggested (^n]) that the 
very large collection area of ground-based imaging at- 
mospheric Cherenkov telescope (lACT) arrays could be 
used to extend CR electron spectrum measurements into 
the TeV domain. The challenge for such instruments (as 
indeed for all CR electron measurements) is to recognize 
electrons against the much more numerous background 
of hadronic CRs. The recent improvements in hadron- 
rejection power achieved by the High Energy Stereoscopic 
System (H.E.S.S.) instrument have now made such a 
measurement possible. 

H.E.S.S. is an array of four imaging atmospheric 
Cherenkov telescopes situated in the Khomas Highland 
of Namibia [H] ■ The array is sensitive to 7-rays (and elec- 
trons) above a threshold of wlOO GeV. The sensitivity of 
the array to extended 7-ray emission has been demon- 
strated with the mapping of supernova remnant shells 
(0j 0)' ^'^'i the diffuse emission around the Galac- 
tic Center [§]. The factor ~10 improvement in 7-ray 
flux sensitivity of H.E.S.S. over previous generation ex- 
periments is based largely on superior rejection of the 
hadronic background. Because this measurement does 
not discriminate between electrons and positrons, elec- 
trons is used generically in the following to refer to both 
particle and anti-particle. The H.E.S.S. electron analy- 
sis presented here is based on the selection of electron- 



like events in regions far from 7-ray sources and sub- 
traction of the remaining hadronic CR background using 
air-shower simulations. The data used were acquired us- 
ing the complete 4-telescope array during 2004 to 2007. 
All data passing quality selection criteria, with zenith an- 
gles smaller than 28°, and targeting extragalactic sources, 
were used in this analysis, amounting to 239 hours of live- 
time. Only the central 3.0° of the ficld-of-vicw was uti- 
lized, with regions within 0.4° of any known or potential 
7-ray source excluded. The energy is reconstructed us- 
ing standard methods. The effective collection area using 
the technique described below is energy dependent and 
reaches «5x 10"* m^ at 1 TeV. The total effective exposure 
of this data set at 1 TeV is therefore «8.5 x 10^ m^ srs. 

The most critical aspect of electron analysis is the ef- 
ficient rejection of the hadronic background. Given the 
relatively high flux of cosmic electrons with respect to 
typical 7-ray sources it is appropriate to make tight selec- 
tion cuts to achieve the best possible signal/background 
ratio. Very hard event selection, including the require- 
ment that all four H.E.S.S. telescopes triggered in the 
event, leads to a greatly increased energy threshold of 
«600 GeV. A Random Forest [9] (see also fioll) approach 
was used to convert image information from the four cam- 
eras into a single parameter C describing the degree to 
which a shower is electron-like. The primary input pa- 
rameters to the Random Forest algorithm are the Hillas 
moments jTl[ of the images recorded in each telescope. A 
C value of zero corresponds to a shower which is almost 
certainly background, and a value of one is assigned if the 
shower is almost certainly an electron. Random Forests 
were trained in five energy bands using simulated electron 
showers and data taken from empty regions. To subtract 
the hadronic background the C distribution of protons 
and nuclei must be known. For this purpose sets of 10^*^ 
proton showers and showers of heavier nuclei were simu- 
lated with CORSIKA ^ using both the SIBYLL ^ 
and QGSJET-II interaction models. About lO'^ 
of these showers trigger the array, and due to the ex- 
tremely efficient background rejection, only 10^^ fall into 
the regime C > 0.9. While a component of heavier nuclei 
is required to explain the distribution of C at values up to 
0.5, the background can be considered as purely protonic 
at larger values of C- 

Fig. [1] shows the measured distribution of the par am- 
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FIG. 1: The measured distribution of the parameter C,, com- 
pared with distributions for simulated protons and electrons, 
for showers with reconstructed energy between 1 and 4 TeV. 
The best fit model combination of electrons and protons is 
shown as a shaded band. The proton simulations use the 
SIBYLL hadronic interaction event generator. The left inset 
shows the complete distribution from zero to one with entries 
on a log scale; the data are shown as points, the filled his- 
togram shows a mixed composition (proton. He, N, Si & Fe) 
cosmic-ray model. To demonstrate the match between sim- 
ulation and data in electromagnetic showers, the right inset 
shows background subtracted 7-ray data as points and 7-ray 
simulations as filled histogram. 



eter ^ compared with the simulated distributions for the 
energy range 1-4 TeV. The peak close to C = 1 is ev- 
idence of a diffuse component of purely-electromagnetic 
showers at these energies. The data at C > 0.6 can be 
described by a combination of simulated electrons and 
protons. By fitting the ^ distribution of the data with 
the distributions of simulated electrons and protons in in- 
dependent energy bands (with two free parameters being 
the electron and proton contribution), the most proba- 
ble number of measured electron showers in each energy 
band can be deduced. The total normahzed goodness- 
of-fit in the C range of 0.6-1 for reconstructed energies 
between 1 and 4 TeV is ~)(^ jv = 0.98 for a model of 
simulated electrons and protons using SIBYLL (proba- 
bility V = 0.5) and 2.15 for a model using QGSJET-II 
(•p = 1.7 X lO^''), which demonstrates that the SIBYLL 
model provides a better description of measurable param- 
eters of air showers initiated by protons of TeV energies. 
Coupled with the knowledge of the energy-dependent ef- 
fective collection area, which is obtained from electron 
simulations following a powerlaw with a spectral index of 
3.3, the number of measured electron showers can be used 
to determine the primary electron spectrum. As electron- 



initiated air showers are in practice extremely difficult to 
separate from 7-ray showers, the peak in our data at 
^ = 1 may contain a contribution from 7-ray showers. 
The signal measured by H.E.S.S. (close to C — 1) is 
therefore a combination of the CR electron flux (CREF) 
and the extragalactic 7-ray background (EGRB). The 
level of the EGRB lies many orders of magnitude below 
the CREF at GeV energies but a naive extrapolation of 
the last few data points measured by EGRET [l^ yields 
an dN/dE oc spectrum which reaches the level of 
the CREF at a few TeV. However, most models for the 
EGRB yield TeV fluxes at least one order of r nag nitude 
lower than this extrapolation (see for example [16|). The 
predicted flux of inverse Compton scattered solar pho- 
tons off CR electrons is also negligible due to our geom- 
etry pointing away from the Sun 17|. Given the uncer- 
tainty in the EGRB/CREF ratio at TeV energies it is 
desirable to separate electrons and 7-rays in our data. 
Essentially the only useful separation parameter is the 
depth of shower maximum (X^ax), which occurs on av- 
erage w20 gcm^^ (or ^ half a radiation length) higher 
in the atmosphere for electrons. Fig. [2] compares recon- 
structed Xmax distributions for simulated protons, elec- 
trons and 7-rays to the experimentally measured ^max 
distribution for electron-like events (C > 0.9). A fit of the 
A^max distribution with the electron/7-ray fraction as a 
free parameter results in a maximum 10% contribution of 
7-rays to the signal (for a confidence level of 90%), which 
is supported by the displacement between the Amax dis- 
tributions from data used for this electron analysis and 
data from a 7-ray rich data set (inset of Fig. How- 
ever, taking into account a conservative systematic un- 
certainty in the determination of Amax of 5 gcm^^ due 
to atmospheric uncertainties, we cannot exclude a sig- 
nificant contamination of ~50% of our electron measure- 
ment by the diffuse extragalactic 7-ray background. Sys- 
tematic uncertainties in the hadronic modeling are not 
considered. 

Fig. [3] shows the CR electron spectrum derived from 
this analysis together with a compilation of earlier 
measurements. Systematic errors on the reconstructed 
spectrum arise from uncertainties in the simulation of 
hadronic interactions and the atmospheric model, as well 
as in the absolute energy scale. The energy scale un- 
certainty is wl5% and is illustrated by a double arrow 
in Fig. [3l The uncertainty arising from the subtrac- 
tion of the hadronic background has been estimated by 
comparison of the spectra obtained using the SIBYLL 
and QGSJET-II models. The C, distributions for pro- 
tons show a slight rise toward C = 1 (see Fig. [T]), pre- 
sumably reflecting events where a large fraction of pro- 
ton energy is transformed to a single tt^ . The rise is 
somewhat more pronounced for SIBYLL as compared to 
QGSJET-II, giving rise to the model dependence. Ar- 
tificially doubling the 7-ray like component in SIBYLL 
reduces the electron flux by ^ 20%, without significant 
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FIG. 2: The distribution of reconstructed shower maximum 
(Xmax) for H.E.S.S. data, compared to simulations. For 
each shower the measured Xmax is corrected for the energy 
dependent shower elongation (93 gcm~^/decade is the re- 
constructed elongation rate expected for electron primaries). 
Showers with reconstructed energies between 1 and 4 TeV are 
included. The bands show the combination of electrons and 
protons (simulated using SIBYLL) and of 7-rays and protons, 
with a ratio determined by a fit to the distribution of the 
data in this energy range. The distributions of electrons and 
7-rays are shown for comparison. The inset contains a com- 
parison of this data (black) with a 7-ray rich data set taken 
from regions < 0.15° from 7-ray sources (gray). 



change in spectral shape. Detailed tests of the analy- 
sis using different zenith angle ranges, different analysis 
cuts (variations of the cuts on the maximum impact 
distance of the showers and the minimal intensity of the 
shower image in the camera) , different regions in the sky, 
different seasons and years as well as another fitting al- 
gorithm all yield consistent results. The estimated sys- 
tematic errors, apart from the 15% scale uncertainty, are 
illustrated by the shaded band in Fig. [31 Our data are 
well described by a power-law: dN/dE = fc(i?/lTeV)~^ 



and 
^ Fit A), 



with k = (1.17 ±0.02) X 10"'' TeV-i m-2 sr-i s" 
r = 3.9 ± 0.1 (stat) [x^/v = 3.6, p = lO^^, 
which implies a steepening of the spectrum compared to 
GeV energies. The spectral index shows little model and 
sample dependence, resulting in Ar(syst.) < 0.3. At 
lower energies the flux reported here is somewhat higher 
than previous results, but fully consistent within the 15% 
scale error. Leaving the scale factor free, H.E.S.S. data 
combined with earlier electron data are well reproduced 
by an exponentially cutoff powerlaw with an index of 
-3.05 ± 0.02 and a cutoff at 2.1 ± 0.3 TeV, combined 
with a scale adjustment of —11% (Fit B). H.E.S.S. data 
are also compatible with very recent ATIC data [2^ , but 



due to the limited energy range no conclusion can be 
drawn concerning the existence of a step in the spectrum 
as claimed by ATIC. 
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FIG. 3: The energy spectrum dN/dE of CR electrons as 
measured by H.E.S.S. in comparison with previous measure- 
ments. The H.E.S.S. data are shown as solid points. The 
two fit functions (A and B) are described in the main text. 
The shaded band indicates the approximate systematic er- 
ror arising from uncertainties in the modeling of hadronic 
interactions and in the atmospheric model. The double ar- 
row indicates the effect of an energy scale shift of 15%, the 
approximate systematic uncertainty on the H.E.S.S. points. 
Previous data are reproduced from: AMS [jal. HEAT [191, 
HEAT 94-95_p3l, BETS [Ml], PPB-BETS [ifTKobayashill 
and ATIC Hj"^ 

Whilst the detailed interpretation of this result is be- 
yond the scope of this paper, we note that our measure- 
ment implies the existence of at least one source of CR 
electrons in the local Galaxy (within ~1 kpc). Some sce- 
narios of a strong local source are excluded. This 
measurement is the first ground-based measurement of 
CR electrons. Future lACT arrays with effective areas 
beyond 10^ m^ should be able to extend the spectrum to 
10 TeV using this technique. 
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